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Abstract 

In the Randall-Sundrum model where the Standard Model fields are confined to the TeV brane 
located at the orbifold point 9 = tt and the gravity peaks at the Planck brane located at 9 = 0, 
the stabilized modulus (radion) field is required to stabilize the size of the fifth spatial dimension. 
It can be produced copiously inside the supernova core due to nucleon-nucleon bremstrahlung, 
electron-positron and plasmon-plasmon annihilations, which then subsequently decays to neutrino- 



antineutrino pair and take away the energy released in SN1987A explosion. Assuming that the 
supernovae cooling rate e < 7.288 x 10~ 27 GeV, we find the lower bound on the radion vev ((f)) ~ 9.0 
TeV, 2.2 TeV and 0.9 TeV corresponding to the radion mass = 5 GeV, 20 GeV and 50 GeV, 
respectively. 
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I. INTRODUCTION 



Several new models based on extra spatial dimensions have been put forward to explain 
the large hierarchy between the Plack scale Mp;(~ 10 19 GeV) and the electro-weak scale 
Mew(~ 100 GeV) Among them the Randall-Sundrum(RS) model is particularly in- 
teresting 2l| since it solves the hierarchy problem in an elegant manner. According to this 
model the world is 5-dimensional and the fifth spatial dimension is characterized by the 
angular coordinate —71 < 9 < 71. The space is an S 1 /Z2 orbifold (i.e. the point (x, 9) 
is identified with the point (2, — 9)). The metric describing such a 5-dimensional world is 
non-factorizable and a line element in this space-time is given by 

ds 2 = e- 2kR ^ v dx»dx v - R 2 c d9 2 (1) 

where k is the bulk curvature constant and x^ are the Lorentz coordinates of four dimensional 
surfaces of constant 9. This theory postulates two D 3 branes along x^ directions living in 
5-dimensional world: one is located at the orbifold point 9 = where gravity peaks(strong) 
and the other at the orbifold point 9 = 71 where the Standard Model(SM) fields reside and 
gravity is weak. The factor e~ 2kRc ^ appearing in the metric is kown as the warp factor. 
The compactification radius -R c (~ the distance between the two D 3 branes) can be related 
to the vacuum expectation value (VEV) of the modulus field T(x) which corresponds to 
the fluctuations of the metric over the background geometry given by R c . Replacing R c by 
T(x), we can rewrite the RS metric at the orbifold point 9 = ix as 

ds 2 = g v ; u s dx»dx u - T{xfd9 2 (2) 

where g™J = e~ 2 ' KkT ^ri ilv = (^j^j V^u- Here f 2 = 24 ^ 5 and M 5 is the 5-dimensional 
^lanck scale. One is thus left with a scalar field <fi{x) which is dubbed as the radion field 
3| . The existence of the radion (modulus) field is a direct and straightforward consequence 



of the existence of the non-factorizable metric. Randall and Sundrum j^j showed that if 
the above metric EqJT] be a solution of the 5-dimensional Einstein equations, then k is to 
be related to the bulk cosmological constant and the vacuum energies of the two D 3 branes 
in a particular way |2j. The modulus field in the original RS model was massless and it 
had no potential: so it was not stabilized. One needs to generate a stable vacuum for T(x) 
at R c , which in turn can give <fi{x) a non-zero vev. This is done in the Goldberger-Wise 
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mechanism [3[, using a bulk scalar field with suitable interactions with the two 3 branes, 
whereupon a potential for the modulus field is generated, and one ends up with a radion of 
nonzero mass. In particular, the radion can be lighter than the other low-lying gravitonic 
degrees of freedom and can very well act as the first messenger of a scenario with compact 
extra dimensions, and reveal itself in collider experiments. Several studies on the observable 
implications of the radion are available in the literature 4]. 

Beside the collider signals, the model predicts a variety of novel signals which can be 
tested in a class of astrophysical or cosmological observations. In particularly, the energy 
loss mechanism of the core-collapse SN1987A explosion and it's relevance in new physics 
context is an exciting area to work. A lot of studies in this direction have already been made 
and are available in the literature |5|. It is interesting to see whether the light stabilized 
brane- world radion do have some role in the supernovae cooling or not. The present work 
is intended to explore this possibility. 

II. SUPERNOVA EXPLOSION AND COOLING 

Supernovae, the final state of an exploding star, come in two main observational varieties: 
Type II are those whose optical spectra exhibit Hydrogen lines and have less sharp peaks 
at maxima (of 1 billion solar luminosities), whereas the optical spectra for the Type I 
supernovae does not have any Hydrogen lines and it exhibits sharp maxima [6]. Physically, 
there are two fundamental types of supernovae, based on what mechanism powers them: the 
thermonuclear supernovae and the core-collapse ones. Only supernovae la are thermonuclear 
type and the rest are formed by core-collapse of a massive star. The core-collapse supernovae 
are the class of explosions which mark the evolutionary end of massive stars (M > 8M ). 
The kinetic energy of the explosion carries about 1% of the liberated gravitational binding 
energy of about 3 x 10 53 ergs and the remaining 99% going into neutrinos. This powerful and 
detectable neutrino burst is the main astro-particle interest of the core-collapse supernovae. 

In the case of SN1987A, about 10 53 ergs of gravitational binding energy was released in 
few seconds and the neutrino fluxes were measured by Kamiokande {t| and 1MB [8] collabo- 
rations. Numerical neutrino light curves can be compared with the SN1987A data where the 
measured energies are found to be "too low". For example, the numerical simulation in 9] 
yields time-integrated values (E Ue ) « 13 MeV, (E Pe ) « 16 MeV, and (E Vx ) « 23 MeV. On 
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the other hand, the data imply (Ep e ) = 7.5 MeV at Kamiokande and 11.1 MeV at 1MB 10 ]. 
Even the 95% confidence range for Kamiokande implies (E Pe ) < 12 MeV. 



lavor oscillations 

would increase the expected energies and thus enhance the discrepancy 10J. It has remained 



unclear if these and other anomalies of the SN1987A neutrino signal should be blamed on 
small-number statistics, or point to a serious problem with the SN models or the detectors, 
or is there a new physics happening in supernovae? 

Since we have these measurements already at our disposal, now if we propose some 
novel channel through which the core of the supernova can lose energy, the luminosity in 
this channel should be low enough to preserve the agreement of neutrino observations with 
theory. That is C newc h anne i < 10 53 ergs s -1 . This idea was earlier used to put the strongest 
experimental upper bounds on the axion mass [11]. In large extra dimension scenario (where 
the weakness of 4 — d gravity is obtained by the large size of the extra spatial dimensions via 
Mp[ = (27T R) d M^~ 2 the KK gravitons interact with the strength of ordinary gravitons 
and thus are not trapped in the supernovae core. During the first few seconds after collapse, 
the core contains neutrons, protons, electrons, neutrinos and thermal photons (plasmons). 
There are a number of processes in which higher-dimensional gravitons can be produced. 
For the conditions that pertain in the core at this time (temperature T ~ 30 — 70 MeV, 
density p ~ (3 — 10) x 10 14 g cm" 3 ), the relevant processes are shown below 

(i) Graviton(^) production in Nucleon-Nucleon Brehmstrahlung: N + N—+N + N + Q 

(ii) Graviton production in photon fusion: 7 + 7 — > Q 

(iii) Graviton production in electron-positron annihilation process: e~ + e + — > Q 

(iv) Graviton production in plasmon-plasmon (photon inside plasma becomes massive and 
called plasmon) annihilation process: 7p + 7p — > Q 

The constraint on luminosity of this process can be converted into a bound on the 4+d 
dimensional Planck scale Mp. Raffelt has proposed a simple analytic criterion based on 
detailed supernova simulations ]5|: if any energy-loss mechanism has an emissivity greater 
than 10 19 ergs g" 1 s" 1 then it will remove sufficient energy from the explosion to invalidate 
the current understanding of Type-II supernovae's neutrino signal. 

The dominant process relevant for the SN1987A where the temperature is comparable 
to and so the strong interaction between N's is unsuppressed. This process can be 
represented as (see above) 

N + N ^N + N + g, (3) 
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where N can be a neutron or a proton and Q is a higher-dimensional graviton. 

Table 1 



Group / Collaboration 


M fl (GeV) 


d 


Cullen et al. [12] 


> 50 TeV, > 4 TeV, > 1 TeV 


2, 3,4 


Barger et al. [13] 


> 51 TeV, > 3.6 TeV 


2,3 


Hannestad e£ a/. [14] 


> 84 TeV, > 7 TeV 


2, 3 



Table 1: The lower bound on the higher dimensional Planck scale Mp> corresponding to the 
no. of extra spatial dimensions d is shown. The lower bound follows from the fact for any new 
physics channel contributing to the SN1987A energy loss, the loss rate e < 7.288 x 10 _27 GeV. 

For the core temperature T = 30 MeV and p = 3 x 10 14 g cm -3 , we list in Table 1 the 
results(lower bound on Md) of various authors. In addition it is worthwhile to mention that 
the KK gravitons produced in plasmon-plasmon collision which contributes in the supernovae 
cooling gives rise strong bound on Mp>: for d = 2 one finds Mp> > 22.5 TeV and for d — 3 
one finds M D > 1.4 TeV [15 1. 

Another extra-dimensional model that can play a crucial role in this SN1987A cooling is 
the Randall-Sundrum(RS) model and was first looked at in 16]. The authors in |l6| studied 
the impact of a light radion on neutrino-anti-neutrino oscillation. They found that for a 
light radion of mass > 1 GeV with ((f)) = 1 TeV, the interaction potential(arising due 
to the exchange of a radion between the supernoave matter and the neutrino-antoneutrino 
pair) does not affect the neutrino oscillation. However, the role of a light radion in the 
Supernovae cooling was not looked at. The present work is an effort in that direction. 

We will see that how the radion produced in electron-positron, plasmon-plasmon annihi- 
lation takes part in the supernovae SN1987A cooling. 




(i) 



k 

->-• 



(ii) 





(iii) 



FIG. 1: Feynman diagrams for the processes 7p + jp — > (f>, e e + — > v x v x and 7p + jp 
v x + ^7 (where x = fx, r) which contribute to the SN1987A cooling. 
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The three primary new physics diagrams that may contribute to the supernovae cooling are: 
(i) 7p + 7p — > 0, (ii) e + + e~ — > (f> — > v x + 77^ and (iii)7p + 7p — > <f) ~ > v x + Vx where x — /x,t. 
These are shown in Fig. [TJ 



III. METHODOLOGY OF CALCULATION 



In Goldberger-Wise mechanism the modulus field T(x) gains a mass and the radion field 
$ ( where $ = Je _7rfcT< - :r - ) with f = Mp;)couples to the brane matter through the trace of 
the energy-momentum tensor [3| 

where <fi = $ — (0) is the fluctuation of the radion field from the VEV ((f)). In the RS 
scenario, ((f)) ~ TeV for kR c ~ 12 in order to produce the weak scale from the Planck scale 
through the exponental warp factor. Since radion coupling to SM matter is determined by 
the 4 — d general covariance, it's coupling to matter is universal i.e. it couples to the trace 
TJf of the energy-momentum tensor T^ v of the SM matter fields which resides on the TeV 
brane and is given by 



T»(SM) = 



Si — — — 



2m 2 w W+W^-m 2 z Z^ 



+ (2m^ 2 - dJWh) + ■■■ (5) 



Inside the supernovae, the relevant matter fields are the nucleons, electrons-positrons, plas- 
mons and the interaction of braneworld radion {(f)) with them is given in Eq. HI 
Now for a generic 2 — > N body scattering, the scattering cross section is given by 



-■LI n ^ -x>) (6) 

where Flux = ^E\E2V re \. Here E±, E% are the energies of the incoming particles 1 and 2 
whose masses are m% and m 2 , respectively and v re i is the relative velocity between them. 
Defining the energy loss per unit mass e = (erg g^ 1 s _1 ) for a generic 2 — > N scattering 
contributing to the cooling process at temperature T one finds 

. 2 r J3„ n r / 2 

(>S X 



^n/^/«)n/^±/^)W' E P .-E ft \ J2 i m i 2 

i=l j =3 \i=l j=3 / spins 

(7) 
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where /j, the occupation numbers for the initial colliding paricles (i.e. electrons, positrons, 
plasmons or nucleons) and (1 ± fj), the Pauli blocking factors for the final state bosons (+ 
sign) and fermions(— sign) 

For a general reaction of the kind a + b — > c, the above expression takes the form 

1 



a 



■\M fi \ 2 27r5(S -m 2 c ). (8) 



Flux 

In the center of mass(c.o.m) frame, we use the notation \/S for the total initial energy 

s/S = E 1 + E 2 (9) 
Flux = 4EiE 2 v rel = 4\p\VS, (10) 

where |p| = |pi| = | p> 2 1 = A 1 ^/^ 1 ' r " 2 ' > an d E x and E 2 are the energies of the particles 1 and 
2. The function X(x, y, z)(= x 2 + y 2 + z 2 — 2xy — 2yz — 2zx), is the standard .fTallen function. 
Now we are concerned with the energy loss to radion. It is convenient and standard 

nn 

[5|, 124J] to define the quantities s a +b^c which is the rate at which energy is lost to radion via 
the process a + b — > c where c has a decay width. In terms of the cross-section a a+ ^ c the 
number densities n a ^ for a,b and the mass density p, e is given by 

(n a n b a( a +b^c)VrelE C m) n 1 s 

£a+b^c. — 1 11 J 

where the brackets indicate the thermal averaging and E cm (= E a + E^) is the center-of- 
mass(c.o.m) energy of the two colliding particles a and b. Note that in the present case, the 
final state radion has a smaller decay width but is stable over the size of the neutron star 



because of it's large life time ~ 10 9 (100 MeV/m) 3 yr (See 23j]) and thus it can escape the 
supernovae while allowing it to cool. 

Next we are interested in the processes in which electron-positron and plasmon-plasmon 
collisions produce radion followed by it's decay to muon and tau neutrino anti-neutrino pair. 
These 2 — > 2 (i.e. a + b — > c — > d + g) scatterings will also contribute in the overall SN1987A 
cooling. The volume energy-loss rate (e) per unit mass of the supernovae matter can be 
calculated as 

1 2 f d 3 r> 4 f d 3 D 1 



2 



12) 
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where /i^, the occupation numbers for the initial colliding electrons, positrons, plasmons 
or nucleons and (1 ± ^3,4), the Pauli blocking factors of the final state bosons (+ sign) and 
fermions(— sign). 



IV. RADION PRODUCTION IN PLASMON FUSION 

Photons are quite abundant in supernovae. Due to plasma effect inside the supernovae, 
photons becomes effectively massive. These massive photons are called plasmons. Our 
interest is in the production of a light stabilized radion <fi in plasmon-plasmon annihilation 

1p(Pi)+1p(P2) ^<P(k). (13) 
The interaction vertex of the plasmon-plasmon-radion 7p(j»i) — r )p(p2) — 4>{k) is given by js]] 

In the c.o.m frame, the 4-momentum vectors associated with the incoming and outgoing 
particles are given by 

p? = (£i,0, 0,p), p£ = (3j,0,0,-p), (15) 
W = (^,0,0,0). (16) 

where p = \p\- It often turns out to be more convenient to keep the polarizations explicitly. 
Here and e° are the transverse and longitudinal polarization vectors of a massive gauge 
boson. For a massive vector boson(e.g. plasmon) with momentum fc M = (E, 0, 0, k) and mass 
ttla, the components of the polarization vector can be written as 

(■+(*) = -^(o,M,o) , (17) 
= ^(0,-1,^,0) , (18) 
<£(*) = —(k,0,0,-E) . (19) 



The polarization vectors satisfy the following normalization and polarization sum conditions 

e^e s ;* = A6 ss \ twW = -V + % e"»efe = 45* . (20) 



m, 

8=1 7 
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The total squared amplitude, averaged over the initial three polarizations(massive plasmons 
have three state of polarizations) states and summed over final states for the process 7p(/ci) + 
7p(fc 2 ) ->■ <p(p) is 



2\ 2 



2+1 



2u 2 



(21) 



Sunstituting this in Eq. El we find the total cross-section cr 7P7P ^ as 



a, 



-\M\ 2 (2n) 



d 3 p 



2(2tt) 3 £ ( 



1 

2~s 



2uV 2 



(0) 



1 - 



2ul 



2-k5 [s — m|) 



(22) 



The volume emissivity of a supernova with a temperature T through this process is obtained 
by thermal-averaging over the Bose-Einstein distribution. Hence, the energy loss rate (e 7p = 
~^Q-y P ) due to plasmon plasmon annihilation is given by (similar to that of the energy loss 
rate via 77 — > vv. 



25|) 



1 



( ^ 7p < -'~7p 7p — > <t> ^* e I E c . m ) 



PSN 
1 1 

PSN ^ J WQ 



e«i/T _ ! 



^ : : cr. 



2wia;2 



7p7p" 



where <7 7p7j ,_>.^ is given in Eq. [221 Note that iV 7p = 4j f^ 1 dui 



1 r oo , ,^(^ 2 -^g) 1/2 • 



,/T_ 1 



is the number density 



of thermal photons or plasmons. In the present case, we treat the plasmon to be transverse 
(with the dispersion relation given by u> 2 = Un + | k| 2 ), since the contribution coming from 
the longitudinal plasmon is typically smaller [26|, |27| . In above ujq corresponds to plasma 
frequency in the supernovae core. 

Finally introducing the dimensionless variables Xj = Ui/T{i = 0, 1, 2) and performing the 
x 2 integration, we find the reaction rate as 



where J 7 



. rxt 1 

^ 9^(<p) 2 p SN 

3 



J-o 



(24) 



+ 



1 / ™±_ 



4 \ %oT 



In above we have used 5 [s — m 2 A 



{x l + x 2 f- r ^ 



(while doing the x 2 integration) where s = {u\ + oj 2 ) 2 = T 2 (x\ + x 2 ) 



Also u>o = Txq where ojq is taken to be equal to (the transverse plasmon mass). 
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V. NEUTRINO PAIR PRODUCTION IN ELECTRON-POSITRON COLLISION 



The neutrino pair production via the s-channel exchange of a light stabilized radion 
(f) produced in electron(e~)-positron(e + ) collision may play a crucial role in the SN1987A 
cooling. The process is 

e~(pi) + e + {p 2 ) A ui(p 3 ) + ui(j> 4 ). (25) 
The radion-fermion-fermion /(pi) — /(P2) — <P(k) is given by [3] 

3« r 8 1 

~W) |#-#"3 m 'J (26) 
where pi, p 2 are the incoming momenta In the c.o.m frame, the momentum vectors for this 
reactions are 

Pi = (Ei, 0, 0, \p\), p£= (E 2 , 0, 0, -\p\), (27) 
p£ = (£ 3 ,0,0,|p*'|), p^ = (E,,0,0,-\p'\). (28) 

The total squared amplitude (averaged over the initial spin states and summed over final 
spin states) for the process is 

Finally, the volume energy-loss rate per unit mass of a supernovae at a temperature T is 
given by 

e = — II / ^f^A/ 2 (l-/3)(l-/4) (27r)V(p 1+ p 2 -p 3 -p4) W? (30) 

which after some rearranging can be written as 

T 7 f f x\x 2 



64ii 5 psn J J [l + exp(x 1 - ^r-)] [l + exp(x 2 + \- 

exp(x\ + x 2 )dxidx 2 



[1 + exp(^ + ^)] [1 + exp(^ _ t^t 



TwrrW ( 31 ) 



where 



m 2 m 2 



T3(0)4 [(T 2 (x! + x 2 ) 2 - m 2 ) 2 + m 2 r 2 ] 
x [T 4 (x! + x 2 ) 4 - 4T 2 (x x + x 2 ) 2 {m\ + mj) + 16771*771* ] 

Here s = (i?i + _E 2 ) 2 = T 2 (xi + x 2 ) 2 and wo is defined above. 

10 



VI. NEUTRINO PAIR PRODUCTION IN PLASMON-PLASMON COLLISION 



The third process that we are interested to look at is the neutrino pair production in 
plasmon-plasmon collision via the s-channel exchange of a light stabilized radion <p i. e. 



lp{ki) + lp{k 2 ) — > vi{pi) + vi(p 2 ). 



(32) 



The total squared amplitude (averaged over the initial polarization states and summed over 
the final spin states) is given by 



M 



iy j2 \m\ 2 = 16m A m »i ( s - 4m ' ; ) 



spins 



2m\ 



(33) 



With this the energy loss rate (volume emissivity) Q due to this process can be written as 



q = II 



2 2 rf 3 k, 



n 



d 3 Pi 



x 



\J 2^(2^11; 2^(2tt) 3 
(2n)W(p 1+ p 2 -p 3 -p 4 )\j2\ M f 



r/i/ 2 (l-/3)(l-/ 4 ) 



(34) 



spins 



where fi i2 , the occupation numbers for the initial state plasmons and (1 — ^3,4), the Pauli 
blocking factors of the final state neutrino and anti-neutrinos. 

Finally, the energy loss rate (volume emissivity) per unit mass at temp. T is given by 
T 7 f f , , x\x 2 



e = 



16n 5 p 



SN 



where 



M 



* / / dx\dx 2 - : : 77 : u _ . -, 

J J [exp{x 1 ) - 1] [1 + exp{x 2 + %-)\ 

exp(xi + x 2 )dx\dx 2 

* [1 + exp(^ + ^)] [1 + exp(^ - ^)] 

16m 2 u m\ [T 2 { Xl +x 2 ) 2 -Ami] 

~ 9T3(0)4 [ (T 2 (a;i+X2) 2_ m 2 )2+m 2 r 2] 



(35) 



2 V 2m\ 



where s = (E 1 + E 2 ) 2 = T 2 (x x + x 2 ) 



VII. NUMERICAL ANALYSIS 



As discussed above the processes of our interest that may contribute in the SN1987A 
cooling are radion mediated three processes: (i) 7p + 7p — >■ 0, (ii) e" + e + v x + and 
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(iii) 7p + 7p — > v x +~Vx (with x — r). If the SN1987A cooling is due to a class of radion 
mediated processes the emissivity rate for those channels must be e < 10 19 erg g~ x s _1 j^J, 
which can be converted to a lower bound on radion vev (0) as discussed below. 



A. Bound on (0) from the + jp 



decay 



In Fig. [2] we have shown the energy-loss rate e{^fp + 7p — > 0) as a function of the radion 
(0) corresponding to an ultra-light stabilized radion of mass = 1.25, 1.5, 1.75, 2.0 and 
2.5 GeV. The upper horizontal line corresponds to the upper bound of the supernovae energy 



1000 2000 30 DO 4000 5000 6000 7000 BOQ0 9000 10000 



FIG. 2: The supernovae energy-loss rate de/dt GeV -1 due to radion emission produced in plasmon- 
plasmon annihilation is shown as a function of {(f)) (GeV). While going from the topmost curve to 
the lowermost curve, increases as 1.25, 1.5, 1.75, 2.0 and 2.5 GeV. The upper horizontal line 
corresponds to the upper bound of the supernovae energy loss rate i.e. e< 7.288 x l(T 27 GeV. 

loss rate e = 7.288 x 10~ 27 GeV for any new physics channel. From the topmost curve to 
the lowermost curve the radion mass increases as 1.25, 1.5, 1.75, 2.0 and 2.5 GeV. We 
find (0) > 7852 GeV for = 1.25 GeV and (0) > 230 GeV for = 1.5 GeV. No bound 
follows for m = 1.75, 2.0 and 2.5 GeV. 



B. Bound on (0) from the e + e + — > v\ + v\ scattering 

The muon and tau neutrino or anti-neutrino produced in e~ e + annihilation (due to a s 
channel exchange of a light stabilized radion) may take away the supernovae energy and 
hence can explain the deficit of electron type of (anti) neutrino at the detector recorded by 
the Kamiokande and 1MB collaborations. On the left panel of Fig. [3] we have shown the 
energy loss rate e(e~ + e + — — >■ + P^) as a function of (0), whereas on the right panel we 
have shown the energy loss rate e(e~ + e + — ^> v T + V^) as a function of (0). The horizontal 
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FIG. 3: The energy-loss rate de/dt GeV -1 is shown as a function of ((ft) (GeV). From the topmost 
to the lowermost curve, rrij, increases as 1, 3, 5, 10 and 50 GeV, respectively. The horizontal 
line corresponds to i < 7.288 x 10 _27 GeV. The left panel corresponds to the energy loss due 

to e~ + e + — > Up + T7j2 scattering, whereas the right panel corresponds to the energy loss by 

, <t> 

e + e" 1 " — > v T + v T scattering. 

line on each panel corresponds to the upper bound e < 7.288 x 10 _27 GeV. On the left(right) 
panel, from the topmost to the lowermost curves increases as 1, 3, 5, 10 and 50 GeV. 



FIG. 4: The lower bound on (<ft) is plotted against m^. The lower and the upper curves follows from 
the fact that the supernovae cooling stems from the e~ + e + v^ + VJl and e~ + e + v T + 
scattering processes, respectively. The region above each curve is allowed. 



In Fig. HI we have shown as a function of ((ft) (the lower bound). The bound follows 
from the fact if any new physics channel contribute to the SN1987A energy loss, it should 
be e < 7.288 x 10~ 27 GeV. The lower(upper) curve follows from the fact that the SN1987A 
cooling is caused by the e~ + e + — > v^Vt) + T^(y^) radion mediated processes. The region 
above each curve is allowed. 

Below in Table 2 we have shown the lower bound on ((ft) for different (refer Fig. 
Hj). From Table 2 we see that for a given channel, the lower bound on ((ft) decreasees with 
the increase in m^. Also for a particular the lower bound on ((ft) (third column) which 
follows from v T — 7v production is much higher than that(second column) follows from v^—TT^ 
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production. Note that the bound on (0) (second column) is much weaker for a radion of 
mass of few tens of a GeV. 

Table 2 



m<£ (GeV) 


(0) (GeV) (lower curve) 


((p) (GeV) (upper curve) 


1 


509 


4846 


3 


166 


1580 


5 


99 


946 


10 


50 


473 


15 


33 


315 


20 


25 


236 


25 


20 


189 



Table 2: The lower bound on (0) (follows from Fig. ^ corresponding to different is shown. 
The second(third) column follows from the energy loss rate e(e~ + e + — — )■ v^t) + ^V(r)) < 
7.288 x 10~ 27 GeV. 

C. Bound on (0) from the + 7p — > v\ + uj scattering 

The Supernovae may cools off by emitting muon and tau neutrino and anti- neutrino 
produced in plasmon-plasmon annihilation inside it's core. The relevant processes are 
7p + 7p — v x + T^x (with x = /i, r). In Fig. [5j we have shown the energy-loss rate 
as a function of (0) corresponding to different values. With the lowering of radion mass, 





FIG. 5: The energy-loss rate de/dt GeV -1 is shown as a function of ((f)) (GeV). From the rightmost 
to the leftmost curve, m^ increases as 1, 3, 5, 10 and 50 GeV, respectively. The horizontal line 
corresponds to the upper bound on the energy loss rate i.e. e < 7.288 x 10" 27 GeV. The left(right) 
panel corresponds to the energy loss by jp + 7p — > Vfi + vjj, (jp + 7P v T + 17^) processes. 
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the lower bound on (0) decreases. In Fig. El we have plotted (0) (the lower bound) as a 
function of radion mass m^. The lower and the upper curves correspond to the SN1987A 



1e+0E 



100000 



5 10 15 20 25 30 35 40 45 50 

rVGeV) 

FIG. 6: The lower bound on (0) is plotted against m<h- The lower bound follows from the fact 
for any new physics channel contributing to the SN1987A energy loss, the loss rate e < 7.288 x 
10~ 27 GeV. The lower(upper) curve follows from the fact that the processes 7p + 7p — > /v M (z/ T ) + 
I^(^v) contribute to the SN1987A cooling. The region above each curve is allowed. 

cooling due to 7p + 7p v^ + V^ and 7p + 7p — ^-f v T + T\ processes, respectively. In Table 
3 we have shown the lower bound on (0) for different values. For an ultra-light radion 
of mass 1 GeV, we find the lower bound on (0) as 44.6 TeV and ~ 425 TeV, respectively 
whereas for a radion of mass about 50 GeV, we find the lower bound as 0.9 TeV and 8.6 
TeV, respectively. 



Table 3 



(GeV) 


(0) (GeV) (left panel) 


(0) (GeV) (right panel) 


1 


44648 


424951 


3 


15042 


143147 


5 


9033 


85958 


10 


4518 


42993 


20 


2259 


21498 


30 


1506 


14332 


50 


904 


8600 



Table 3: The lower bound on (0) as a function of corresponding to the SN1987A energy 
loss rate e < 7.288 x 10~ 27 GeV. The second(third) column of the lower bound on (0) follows 
from the energy loss via the channel 7p + 7p — v^t) + ^V(t) ■ 
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VIII. CONCLUSION 



In this work we have investigated the impact of a light stabilized radion in super- 
novae cooling. The radion produced inside the core of a supernova due to electron- 
positron and plasmon-plasmon collisions, can take away much of the energy released in 
supernovae explosion. The primary processes of our concern are (i) 7p + 7p — > 0, (ii) 
e + + e~ z/ M(r ) + I7 M ( r )and (iii) 7 P + 7 P v^ T ) + V^y Assuming that the energy 
loss rate due to each of the above three channels e < 7.288 x lCT 27 GeV, we obtain the 
following lower bound on the radion vev (0) due to a light stabilized radion. For the pro- 
cess 7p + 7p ->• 0: we find (0) > 7.85 TeV for m = 1.25 GeV and (0) > 2.3 GeV for 
= 1.5 GeV. For e+ + e~ -A u x + V x : with x = t we find (0) > 0.95 TeV and > 0.24 
TeV corresponding to = 5 GeV and 20 GeV, respectively. No reasonable bound on (0) 
follows with v^Vp as final state particles. Finally the process 7p + 7p — v x + V x : for 
x = /z, we find (0) > 9.0 TeV and > 2.4 TeV for = 5 GeV and 20 GeV, respectively. 
For x = r, we find (0) > 86 TeV and > 21.5 TeV corresponding to the same values. 
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